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Abstract 


The contont of this Quarterly Report can be summarized as follows: 

a) The tether control law to retrieve the satellite has been modified in 
order to have a Bjnooth retrieval trajectory of the satellite that minimizes the 
thruster activation. 

b) The satellite thrusters have been added to our rotational dynamics . 
computer code and a preliminary control logic has been implemented to simulate 
them during the retrieval maneuver. 

c) The high resolution computer code for modelling the three-dimensional 
dynamics of untensioned tether, SLACK3, has been made fully operative and a set 
of computer simulations of possible tether breakages has been run. 

d) The distribution of the electric field around an electrodynemic tether 
in vacuo severed at some length from the Shuttle has been computed with a 
three-dimensional electrodynamic computer code. 
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1.0 INTRODUCTION 


This is the third quarterly report submitted by SAO under contract 
NAS8-36160, "The Investigation of Tethered Satellite System Dynamics," Dr. 
Enrico Lorenzini , PI , and covers the period from 15 February 1985 through 14 May 
1985. 


2.0 TECHNICAL ACTIVITY DURING REPORTING PERIOD AND PROGRAM STATUS 


2.1 Retrieval Control Law 


The existing version of the DUMBEL program including the rotation of the 
subsatellite does not have a retrieval control law, but is easily modified to 
include various control laws. For this project a modified version of the 
control law used in program TETHER has been incorporated into DUMBEL. The 
TETHER program developed at Marshall Space Flight Center uses a tension control 
lav; of the form: 

r - m.(w2(f-f 0 ) + (l + BOOM) (2.1.1) 

+ 2£ c We (4“^ 0 ) ) 

where 


r - tethf?r tension 

m. - effective mass (tjisat + TOteth/2) 

£ c - control law damping 

w 0 = orbital frequency 

BOOM = boom length 

l — tether length 

6 — B ~.0003StQ 

c c = t 0 e + BOOM 

l a = - , 00035l c 


( 2 . 1 . 2 ) 


neglecting the effect of the boom, the control frequency is 
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Wo « 2 w 0 (2.1.3) 

The time t 0 Is the time measured from the start of retrieval. The first 
term keeps the retrieval on a specified time schedule, the second provides the 
equilibrium tension required to support the subsatellite, and the third provides 
damping and controls the rate of retrieval. The control law is used in 
conjunction with thrusters to allow rapid retrieval. The law may be used with 
arbitrary initial conditions. For a 20 km tether, the initial commanded 
retrieval rate is 700 cm/sec. 


2.1.1 Retrieval At A Constant In-Plane Deflection Angle - 

The equation of motion for the in-plane deflection, if we neglect the 
tether mass, is 

Eo — mj ($£ cos ^ 2 ($ + w„) (£ cos sin $) -f 3£wJ cos 0 cos $ sin 0) (2.1,4) 

Where mj is the satellite mass, If we have <f>=>^-d~Eo-0, 
then the equation of motion reduces to 

0 *— mj(2 u 0 i + 3tu} cos 6 sin 6) 


Solving for the retrieval rate we have 


i = — ~ lw„ cos 6 sin 0 


i — — at 


where 


a — cos 6 sin 0 


(2.1.5) 

( 2 . 1 . 6 ) 


(2.1.7) 


The distance £ as a function of time is, by integration 
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t « £„ o' ot (2,1.8) 

Tho acceleration is 

l — al ** a 3 l (2.1,9) 

The equation of motion of the radial variable is 

r = mj {i - ifr - t cos 2 ^ (^ + Wo) 1 + tup o - 3£u»2 cos 2 0 coa^) (2.1.10) 

I f we have <f> = ^ « 0 = 0 , the equation reduces to 

r = m 3 (/ - 3£«>2cos J 0) (2.1.11) 

In order to implement a tension control law we need to evaluate l in order to 
calculate the tension r. Eor retrieval at a constant angle 6, the acceleration 
is, from equations (2.1.9) and (2.1.7) 

i s= a 3 l =s cos 2 0 sin 2 0£ (2.1.12) 

Putting equation (2.1.12) into equation (2.1,11) gives 

T = raj (—wj cob 2 0 sin 2 0£ - 3£w§cos 2 0) = 3m 3 £w§ cos 2 0 (^ sln 2 0 - 1) (2.1,13) 

4 4 

In the equation for tho tension, the distance £ should be measured from the 
orbital center of the system, For a system of masses mj, = Shuttle and m 3 = 
satellite, the distance of m 3 from the center of mass of the system is 
£mi/ (mi + m 3 ) . 

Eor a simple harmonic oscillator of frequency w, the critical damping 
coefficient is 

b = 2mw (2.1.14) 

Eor the in-plane variable 0, the frequency is \/3u 0 . 


Eor the out-of-plane 
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variable 4 the frequency is 2o> 0 and the critical damping coefficient is 4mw„. 

In equation (2.1.4) we sao from the last term that the restoring term is 
approximately 3m£w20 for small angles. The stiffness is therefore 3mw2, since 
tO is the in-plane displacement , For the out-of-plano variable the stiffness 
obtained for small displacements is . 

In the control law of equation (2.1.1) the first term is designed to keep 
the motion along the trajectory proscribed by the commanded length. The 
stiffness is the same as the out-of-plane restoring stiffness. The last term is 
designed to provide damping and uses the coefficient for critical damping of 
out-of-plane oscillations. The middle term is designed to provide the 
equilibrium tension necessary to balance the gravity gradient. 

The retrieval profile given in equation (2.1.2) will produce an in-plane 
displacement which can be calculated from equation (2.1.7) . With o> 0 = 
.00113137, the retrieval angle for a = .00035 is 12. IS degrees. The initial 
retrieval rate for a 20 km tether is al 0 = 7 m/sec. 

The first term in equation (2.1.1) is not necessary for maintaining the 
stability of the retrieval. As an alternative to forcing the retrieval to 
"catch up" to a desired time table, one could omit the first term and allow the 
retrieval to proceed in a natural way and determine the time required for 
retrieval from a hanging position. The difference between this time and the 
time calculated from the commanded length can *ce used as a correction factor to 
adjust the start time of the retrieval to t'ijtain a desired termination time. 
The first term can have the effect of forcing the retrieval in a way that may 
not have a beneficial effect on stability. The value of i c in the last terra 
does not necessarily have to be computed from t c . It could be computed as 
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l a ss — at (2 . 1 • 15) 

This should promote stability by bringing the retrieval reto to the value 
necessary for maintaining the desired retrieval angle for the actual length of 
the wire at any particular time. If the rft&rieval is on a perfect trajectory 
the first and third terms will bo zero since t ~ i a and i ® i a . The second term 
should then provide the tension necessary to maintain a perfect retrieval 
trajectory. The second term in equation (2.1.1) is obtained by sotting l - 0 ~ 
0 in equation (2.1.11). A more accurato expression can bo obtained by taking 
into account the deceleration, the dependence of the tension on the in-plane 
angle, and the correction for the location of the center of mass of the system. 
Th'c> center of mass is not exactly the same as the point of zero radial 
acceleration but the difference can be neglected for this purpose, An 
alternative control law taking into account the above considerations is 

r = 3m«fw§ cos 1 0 c ('j sin 3 (? 0 - l)mi/(mx + mj) + 4w„(i-i c ) (2.1.16) 

where 

L = -al 

3 3 

a — t-w,, cos0 c sin<? c = ■TW o 6in(20 o ) 

2 *1 

In equation (2.1.16) the angle 6 e is a fixed parameter of the tether control 
law. The value of 0 O determines the duration of the maneuver. It also 
represents the Inclination of the steady state trajectory followed by the 
satellite during retrieval. 
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2,1,2 Software Modi ficatitms - 

There are two versions of DUMBEL. One has rotation of the subsntollite and 
the other does not, The version with rotation does not havo many of the special 
features included in the other vorsion. In order to have an initial in -plane 
displacement in the rotational version it was necessary to add soma subroutines 
from the other version. This turned out to be more difficult than anticipated 
because the numbering of the masses in the rotational version is reversed from 
that used in the other version and in SKYHOOK. For consistency it was decided 
to change the numbering systom in the rotational vorsion to agree with that used 
in the other programs, Changes had to be made in programs COORD, D1ER0T, 
DUMBEL, ROTCOORD, and SETUP. Subroutines INITIAL, EQUIL, and LAUNCH have boon 
added to the rotational version so that runs can be done with initial 
displacement in the in-plane and out-of-plana directions, and initial retrieval 
velocities . 

In order to be able to plot the orientation of the subsatellite in various 
coordinate systems, the output of the program has been changed to include the 
rotational part of the state vector. The postprocessor haB been modified to 
plot the orientation of the subsatellite in the rotating orbital coordinate 
system. This was done by adding subroutines ROTCOORD, ELF.M, and ANGROT to the 
postprocessor which has been renamed RSTAVEC. 

A second postprocessor PLOTFIL has also been modified to provide additional 
information. Program PLOTFIL plots the in-plane, out-of-plane and radial 
components of the position of the satellite. For constant wire length, the 
ln-plane component gives a good representation of the in-plane angle, However, 
during retrieval the angles must be computed at each output point from the 
components since the wire length is constantly changing. The program has been 


Pago 10 


changed to plot tho in-plane <n*4 out-of-plano angles ao that the retrieval can 
be monitorod to see If tho retrieval angle la constant. 

Implementation of a tansion control retrieval algorithm in DUMBEL ia 
relatively otraight forward. The tension is computed in subroutine SETUP, Tho 
control laws given by equation (2.1,1) and equation (2.1.16) have boon 
implemented with control parameters tho.t allow various combinations of terms to 
be used, Tho program can be run in either the retrieval mode or tho steady 
state mode. Various debugging runs wore done until tho program appeared to bo 
giving correct results , One bug in particular kept making the run go unstable 
very quickly by applying too much tension. 


2.1.3 Computer Simulations - 

A set of computer simulations has been done to test the retrieval control 
algorithm in DUMBEL. The software modifications described in the previous 
section have been done concurrently with the test runs. Tho following 
parameters have been used in the runs, A 550 kg 3ubsatellite is deployed upward 
on a 20 km tether from a 100 metric ton Shuttle in a circular orbit at 400 km. 
(Atmospheric drag is not included in these simulations, so that the results are 
not significantly dependent on altitude,) 

An initial test run was done with rotation but the run was not completed 
because the integration was too slow. Using the actual moments of inertia of 
the subsatellite gives a relatively short period for rotation under the 
restoring torque of the tension in the wire. This slows down the numerical 
integration. The run was terminated at about 4000 seconds of orbital time with 
the subsatellite at 4.8 km from the Shuttle with an in-plane displacement of 
12.42 degrees. Equation 2,1.16 was used without tho correction for center of 
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mass . 


Since the rotational dynamics is not necessary for studying tho ratrioval 
algorithms it was decided to put the rotrieval facility into tho version without 
rotational dynamics which has the advantage of rapid numerical integration. Tho 
first run with tho program was done with tho control law of equation (2,1,16) 
without tho center of mass correction, Tho run was terminated at 14,000 aoconds 
of orbital time with tho aubsatellita 100,9 meters from the Shuttle, The 
distance oxpocted in a perfect rotrievol at this time UBing equation (2,1,0) is 
112 motors, Tho in-plano angle at tho end of tho run was 12.44 dogreos , Tho 
Initial conditions for tho run were with tho tether diBplacod 12,1046 degroos in 
tho in-plane direction and a retrieval velocity of 7 metors/soc. Plots of tho 
liberation angles showed no out-of-plane displacement, and a decaying in-plano 
oscillation converging to 12,44 degrees. figure 2,1,1 shows a plot of tho 
In-plan,?, ang?,® in degroos vs. time in seconds. 

The second test run was done with tho rotational version of DUMBEL. Since 
tho slow integration is caused by the short poriod rotation of the subsatollite, 
the moments of inertia of the subsatellito were arbitrarily increased by o 
factor of 1000 in order to see if rapid integration could bo achieved for tho 
long runs required during retrieval. Tho run was terminated at lc3,200 second 
with the subsatellite at 32 meters. The integration proceeded rapidly. The 
distance from the Shuttle expected from equation (2.1,8) is 34 meters. Tho 
final value of the in-plane angle was 11.59 dogreos. A plot of tho ln-plano 
angle shows oscillations which are not damped and appear to be increasing at tho 
end of the run. Figure 2.1.2 shows a plot of the in-plano angle in degrees vs. 
time in seconds. The algorithm does not at present treat tho difference between 
thfi attachment point and the center of mass of the satellite in a completely 
rigorous manner. The rotational dynamics may be affecting tho retrieval 
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slightly. 

The third test ran was done using the control law of equation (2.1,1) with 
BOOM = 50 meters, and = 1. The run went unstable and the in-plane angle 
reached 90° at 8500 seconds, The wire length was 1,6 km. The run was done 
without rotational dynamics. The initial conditions have the tether at 12,18 
degrees from the vertical in the in-plane direction and an initial retrieval 
velocity of 7 meters/sec. 

In order to test the effectiveness of the damping in the algorithm for 
removing initial transients , a fourth run was done starting with the wire 
vertical and the subsatellite stationary. The run was terminated at 21,600 
seconds with the subsatellite at 11.32 meters (since the boom length = 12 meters 
this simulates a complete retrieval maneuver) . The value calculated from 
equation (2.1.8) is 10.4 meters, The oscillations of the ln-plane angle were 
well damped, decreasing by about a factor of 2 on each oscillation. The final 
angle was 12.44 degrees. Figure 2.1.3 shows a plot of the in-plane angle in 
degrees vs. time in seconds. The run was without rotational dynamics. This 
run was very well behaved in contrast to the third run. However the final value 
of the in-plane angle is clearly different from the value of 12.1846 used in 
computing the control law constants. The run used equation (2.1.16) without the 
center of mass correction. 

For the fifth run, the correction factor for the center of gravity of the 
system was included in equation (2.1.16). The tether was Initially displaced 
12.18 degrees in the in-plane direction and the initial retrieval velocity was 7 
m/sec, The run was terminated at 21,600 seconds with the subsatellite 10.28 
meters from the Shuttle. The theoretically expected final distance is 10.42 
meters. The small oscillations present initially in the in-plane angle were 
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well damped and the final in-plane angle was 12.1043 degrees, in close agreement 
with the theoretically expected value of 12.1046. 

Three more teat runs have been done in an effort to identify the causa of 
the instability in the third run. In the sixth run, equation (2.1.16) was used 
along with the first term of equation (2.1.1). The boom length we i set to zero. 
Initial conditions were with the tether displaced 12,1046 degrees in-plane and 
an initial retrieval velocity of 7 maters/sec. The run was terminated at 10,700 
seconds with a final wire length of 463 meters. The theoretical value of this 
time is 472 meters. The in-plane angle was stable and the final value was 12.26 
degrees . 

The seventh run was done with the second and third terms of equation 
(2.1,1) and the same initial conditions. The run was terminated at 10,100 
seconds and the final wire length was 178 meters. The theoretical value is 583 
meters at this time. The final angle was 17.16 degrees, It is apparent from 
the results that too much tension is being applied by the algorithm and the 
result is a too rapid retrieval with a large In-plane displacement. 

The eighth run was done with the same parameters as the third run but 
setting the boom length equal to zero. The run was unstable at the same point 
as in the third run. 
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2.1.4 Conclusions - 

The results of the computer simulations show that it is possible 
(neglecting effects such as atmospheric drag) to doviso a control law using only 
length, and length rate information that will bring in the subsatellite along a 
straight line at a desired angle with good damping of in-plane oscillations 
(without any in-plane thrusters) . The success of the control law depends on 
accurately modelling the tension necessary to keep the system in equilibrium at 
the desired retrieval angle, Effects that must be accounted for are the 
dependence of the gravity gradient on the in-plane angle, the deceleration of 
the subsatellite, and the position of tba center of mass of the system. The 
control law has little effect against out-of-plane oscillations which can be 
handled by thrusters. The effects of atmospheric drag can also be handled by 
thrusters . 


2.2 Implementation Of Thrusters In Rotational Dynamics Program 

The TSS currently under consideration will have thrusters in the ln-plane, 
out-of-plane, and In-line directions, and another thruster for controlling 
rotation about the vertical axis of the satellite, It is assumed that the 
tension in the wire will align the subsatelllta attitude along the direction of 
the wire. The thrusters will be used to control llbrations during retrieval. 
In order to Implement the thruster algorithms it is necessary to construct unit 
vectors in the direction of firing of the thrusters and calculate the swing 
velocities in the in-plane and out-of-plane directions. 
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2.2,1 In-Plano And Out-Of-Plano Unit Vectors - 

There are a number of coordinate systems involvod in simulating the use of 
thrusters on-board the subsatellito , The position, velocity, and orientation of 
the subsatellite, and the position and velocity of the Shuttle are integrated in 
inertial coordinates. Since the thrusters are attached to the subsatellite, the 
direction of firing is specified by the direction cosine matrix giving the 
subsatellite orientation. For this initial analysis it is assumed that the 
orientation of the subsatellite is controlled so that the thrusters point in the 
in-plane, out-of-plane, and in-line directions. 

The first task is to construct unit vectors in the direction of the 

thrusters. If Pi and Vi are the position and velocity cf the Shuttle, the 
» 

normal to the orbit is 

N = Pi x Vi (2.2.1) 


The direction of the wire is 

w o P a - Pi 

where P a is the position of the subsatellite. The direction of the in-plane 
thruster is given by 

UiH s= N x w. 


The out-of-plane vector is 


Uoux — w x Gin 
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2,2,2 Libration Velocity Of The SubBatollite - 

The objective ot the use of thrusters during retrieval is to eliminate 
librations that could build up and cause unstable behavior. The input to the 
control algorithm is the libration velocity, In order to calculate this 
quantity it is necessary to know what the velocity of the subsatellite would be 
in inertial space if the subsatellite were not librating. The difference 
between this and the actual velocity is then the libration. The components of 
this velocity along the in-plane and out-of-plane vectors give the in-plane and 
out-of-plane libration velocities. 

If there is no libration, the whole tethered system rotates like a rigid 
body with the normal to the orbit as the axis of rotation. The unit vector 
normal to the orbit is 

ft = N/| N | 

where N is given by equation (2.2.1), The velocity of the subsatellite is in 
the direction 

Q = n x Pj 

The radius of rotation of the subsatellite is 

pi = IqI 

The radius of rotation of the Shuttle is 

Pi - I Pi I 


vi = | Vj. | . 


and the velocity of the Shuttle is 
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The ratio of the subsotollite volocity to the Shuttle velocity is pj/pt without 
libration. For a circular orbit the velocity of the subsatellite should be in 
the direction of Q, However, in an eccentric orbit, there will be a radial 
component of the velocity except at perigee and apogee. In an eccentric orbit 
there is no uniformly rotating orbital reference frame so that it is not clear 
what the velocity of the subsatellite should be without libration. However, in 
an attempt to define an approximate rotating reference system, the radial 
velocity of the Shuttle has been accounted for in calculating the direction of 
the subsatellite velocity without libration. The angle by which the Shuttle 
velocity deviates from that of a perfect circular orbit is approximately 

a ~ Pi.Vi/(/3iVi) 

The unit vector for the direction of the subsatellite velocity is then 
constructed as 

$ = Q /p2 + atPi/pi 

The magnitude of the velocity of the subsatellite is 

v 3 » ViPa/pi 

so that the vector velocity for the subsatellite without libration is 

v a = v 3 q. 

The libration velocity is then 

AV 3 = V 3 - v 3 

where V 3 is the actual velocity of the subsatellite. The in -plane component of 


the libration is 
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V,„ « AVa'Uju 


and the out-of-plane component is 

Vout • Bt AVj'Uout 

There will of course also bo a radial component of AV but this is not needed for 
the thruster control algorithm. 

2.2.3 Software Implementation Of Coordinate Systems For Thruster Algorithms - 

The equations developed in sections 2.2,1 and 2.2.2 have been coded in 
subroutine INOUT (Z, DV, T, WIRE, UIN, UOUT) , when Z is the state vector for the 
Shuttle and subsatellite , DV is the libration velocity in inertial coordinates, 
T is the time, and WIRE, UIN and UOUT, are unit vectors in the radial, in-plane, 
and out-of-plane directions. 

A calling program has been written to read Z and call INOUT for testing 
purposes. A debugging version of the program was written to print the various 
vectors computed. Tests have been done with no libration angle, an in -plane and 
an out-of-plane displacement of the subsatellite for a circular orbit. All 

quantities have been checked against hand calculations. The tests have been 
repeated with a slight eccentricity to see that the vectors change in an 
appropriate manner . 

Since the tests of the subroutine seemed satisfactory the subroutine was 
put into the DUMBEL program. It is called from subroutine DIEFUN which computes 
the time derivative of each component of the state vector of the system. A 
number of test runs have been done and the debugging print examined to make sure 
the results are reasonable. The first test run was with no libration. The 
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differential velocity was zero no it ehould bo. Huns wora dono with on in-plane 
displacement only, and an out-of-plana displacement only, in an equatorial 
orbit. The DUMBED progrom has a simplified facility for sotting up initial 
libration nnglos written for equatorial orbits only. A toot run was dono at 28° 
inclination with no libration, Tho differential velocity was zero as it should 
bo. A final run waa done at a 45° orbital Inclination with on ln-plano 
displacement introduced manually in the initial conditions, The output shows an 
in-plane component of tho differential velocity, no out-of-plane component and a 
slight radial component. 


2,2,4 Critical Damping Of Librations With Thrusters - 

For a simple harmonic oscillator the critical damping coefficient is 2mw, 
where w is the frequency and m is the mass, For in-plane librations, the 
frequency of oscillation is whore w 0 is the orbital frequency. For 

out-of-plane oscillations w = 2w 0 * The critical damping coefficients are 

therefore 

biw — 2mj\/3 w o 

and 


b 0 ui 2mj2tJo — 4m3Wo , 

The thrust required to provide critical damping Is, for each component, 

Fin “ - b iN V IN 


Four — — b 0 ui Vouj. 
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2,2,5 Computer Simulations Of Thrusters for Critical Damping - 

A couple of tost runs have been done using the algorithms for critical 
damping. In tha first run tho tether was given an out-of-plano displacement of 
10 degrees, Tho simulation was run for 5000 seconds, The out-of-plano 
displacement approached tho rest position exponentially without overshooting, 
Tho times corresponding to 10°, 1°, .1° and .01° wore 0, 1700, 2900, 4100 

seconds, Tho amplitude decays by a factor of ton every 1200 seconds, Because 
of the coupling between the in-pleno and out-of-plano angles, tho out-of-plano 
motion Induced an in-plane oscillation of about ±,B°, This in-plane oscillation 
was not damped. Fig, 2.2.4 shows tho out-of-plano angle vs, timo, and Fig, 

2.2.5 shows the in-plane angle vs. time. 

In tho second test run, tho tether was given an in-plane displacement of 
10°, There was no out-of-plane oscillation as expected, since the in-plane 
oscillation does not couple to the out-of-plane as long as the out-of-plano 
displacement is zero. The in-plane angla returned exponentially to zero, Tho 
times when the angle was 10, 1, .1 and ,01 degrees were 0, 200, 3400, and 4800 
seconds. The time to denay by a factor of 10 is about 1400 seconds. Figure 

2.2.6 shows the in-plane angle vs. 


time , 
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Figure 2.2.5 In plane angle (degrees) vs, time (sec.) starting with n 10° 

out-of-plane displacement and using the out-of-plane thrusters for critical 
damping. 
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Figure 2,2.6 In-plane angle (degrees) vs. time (sec.) starting with a 10 c 
in-plane displacement and using the in-plane thrusters for critical damping. 
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L.3 Slack Tether Studios 

The slack tether simulation program SLACK3 has been enhanced to make the 
three dimensional version fully usable, and to improve handling of difficult 
cases. Several case studios have been run and are presented here. 

The influence of tother viscous damping is examined and found to be 
possibly significant, both for the physical tether and our attempts to model it. 


2.3.1 Slack3 Modification And Utilization - 

During the reporting period SAO has enhanced and modified the SLACK3 slack 
tether simulation program, A series of case studies suggested by Mr. C. C. Rupp 
of NASA ware run. Difficulties encountered with these runs led to further 
modifications to the bounce handling routines. 

2. 3. 1.1 Program Modifications And Enhancements - 

As features (such as boom rotation) were added to SLACK3, the number of 
input parameters increased also, until the arcane input procedure required 
approached in complexity that of SKYHOOK, In particular, interactive runs 
involved a long dialogue with much room for operator error and frustration. In 
the current version of SLACK3 , a set of lexical routines developed for SKYHOOK 
(under NASA Grants NA.G5-325 and NAGS-458) were used to create a "menu" 5 most 
parameters now have default values (e.g. a 1/2 tonne subsatellite, deployed, 
upward 20 km) which are displayed Initially; if the user wishes to change, say, 
the angle at which the tether is deployed she will then enter and identifier 
"tether . angle" ; if the identifier is followed by the appropriate number of 
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parameters (e,g. ! 'tathor . angle -35,0 0.0") they are accepted, and otherwise the 
program prompts for the specific quantities. In addition to being "user 
friendly" for interactive use, this structure makes the input file for batch 
processing much more readable, since each non-default parameter is clearly 
identified. Moot input is also done in appropriate units, o.g. km instead of 
cm, degrees Instead of radians. The tether properties can be entered either as 
simple material properties (such as E) or as specific tether properties (AE) . 

Other enhancements to SLACK3 are : 

- The boom origin need no longer be at the Shuttle center of gravity. An 
(x,y,z) may be specified. 

- Thrusters on board the Shuttle may be fired. Rather than "operating" the 
Shuttle's particular thruster configuration, the user must specify an 
acceleration value, direction and initiation time. 

In addition to the enhanced features, several modifications were made to 
SLACKS to avoid "infinite loop" type situations encountered while running the 
case studies described below. Some of these studies involved, for instance, 
long time scales for which the program had not been tested and which led to 
problems with the machine roundoff. The root finding routine, REEINE, was made 
substantially more robust: the interval error tolerance is checked to ensure 

that it is not smaller than the machine's precision (a problem possible to run 
into when using an absolute input tolerance and the time value becomes larger 
than previously experienced) ; corrective action is taken if the Interval does 
not shrink; and as a last resort, the problem .Is solved with the inefficient 
but infallible bisection method (the major algorithm is modified regular falsi ) . 
The method of choosing the left (^arly) limit of the interval containing the 
root in BOOMTIME was analyzed and modified to avoid the occasional case where an 
interval not containing the root was chosen (due to roundoff, root- finder 
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truncation, etc.) whan finding the next bounce between the boom and mass 1 after 
a given such bounce; the default left limit in such a case is the actual bounce 
time, and small errors make this time not specifiable precisely so that one must 
check and occasionally adjust this left limit. One case study, which as 
discussed below is actually physically unsuitable for SLAC5C3 simulation, still 
leads to an infinite loop type of failure; though this is a symptom of the 
underlying physics, we shall attempt to improve the program so a3 to avoid 
outright failure, 


2. 3. 1.2 Case Studies Of Tether Break Situations - 

At the suggestion of Mr. C. C. Rupp of NASA, we have used SLACK3 to study 
several specific examples of operational Interest. All cases involved an 
electrodynamic tether deployed straight up, with the boom deployed 35° backward 
from vertical. No avoidance maneuvers {rotation, thrusters) were made. Tether 
properties of d = 0.3 cm, (i = 0.08 g/cm, and AE = 0.6 x 10 10 dynes were used; 
the damping was set to zero since the correct value is uncertain and damping may 
have serious ramifications for the SLACK3 model (see Section 2.3.2 below). The 
initial configurations were slightly perturbed, the same relative perturbations 
being employed in each case. Each case was run for approximately the same 
physical time relative to a crude time scale defined by the cut length divided 
by the recoil velocity. 

Two series of cases were simulated: Eirst, a set with the tether fully 
deployed to 20 km and severed at 1, 10 and 20 km. Second, breaks at 0.5 km from 
the Shuttle with original deployed tether lengths of 1 km, 10 km and 20 km. For 
the 1 km original tether, a 2 Newton thruster was included by using a tension 
"fudge factor" of 2.0 (since the gravity gradient tension is about 2.08 N) 
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doubling the recoil velocity. 

The results of these simulations are shown in Figures 2.3.1 through 2.3.5, 
in the familiar SKYHOOK side-view and front-view series of snapshots. 

In Figure 2.3.1, the 20/1 case (20 km deployed tether cut at 1 km from the 
Shuttle) , the recoil velocity is relatively rapid, and we see that initially the 
tether is brought forward of the Shuttle by Coriolis force, largely missing the 
Shuttle) after passing the Shuttle and becoming extended in the downward 
direction, air drag forces and an overall "bounce" from the initial forward 
direction take over and the configuration becomes trailing, It is not clear 
whether this trailing position is permanent or whether there will be another 
overall bounce bringing the tether back forward to impinge on the Shuttle; one 
case described in the first quarterly report indicated that for at least some 
tether lengths the trailing position is stable. If of Interest, this case could 
be carried further. 

The behavior shown in Figure 2.3.2, the 20/10 case, is surprisingly 
different. The tether crumples up slightly on recoil, then gravity gradient 
forces take effect and straighten it out again only to have the elastic tether 
"bounce" back in a second recoil, and so forth. The simulation did not go 
nearly as far relative to the recoil time scale as that in Figure 2.3,1 because 
the attempts to come into tension at the overall "bounces" greatly increase the 
number of ball-and-spring bounces which in turn determines the computational 
effort) this case terminated on exceeding a cpu time limit. The Shuttle is 
clearly safe from the vast majority of the tether In this case, since most of it 
does not have the energy to overcome the gravity gradient force, but it could 
become entangled with the portion near itself. This might be clarified by a 
simulation with segments strongly bunched in the region near the Shuttle, and a 


RflOIBL COMPONENT IWETERSl 

6Q. OO -QO.DO -?0-C3 -0.00 





RHQIRL COMMENT IMFIERS1 t *lO ^ 



Pago 33 


modification of the plot routine to allow plotting only the first few segments. 

There is no figure for the 20/20 case. The bouncing behavior due to 
insufficient recoil velocity noticed in the 20/10 case would be, of course, oven 
more pronounced, and led to an Infinite bounce cycle. This failure is being 
studied, but the overall behavior should be about as in the 20/10 case. 

Figure 2,3.3 shows the 20/, 5 case. This is similar to the 20/1 case, with 
the same recoil velocity but even less tether remnant. The overall forward 
contribution of the Coriolis force is not so pronounced (since it does not have 
as much time in which to act) , hence neither is the rearward bounce; the drag 
also does not have as much time in which to act. Figure 2.3,4, the 10/. 5 case, 
is similar, with less forward Coriolis (hence less bounce back) and more drag, 
due to the decreased velocity and lengthened time scale. Interaction * with the 
Shuttle is not clear in either case but seems likely. 

Figure 2.3.5, the 1/.5 case, shows the result of a comparatively slow 
recoil: drag dominates completely, At least for the initial recoil the tether 
stays as clear of th:j Shuttle as one could desire, though there is some 
indication that the tether may bounce forward, leading to impact. This case 
terminated due co excess cpu time: the dominant drag tends to bring the tether 
into tension, leading to frequent ball-and-spring bounces. 








Figure 2,3.4. 10 km tether cut at 0.S km. Output at 25 second intervals, total 

run of 1450 sec. 
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2. 3, 1,3 Pallui'y** And Difficult Cases! Implications For Tha Model - 

Those cases which either foiled or took inordinate amounts of computer time 
due to largo numbers of bounces, both in the above simulations and from provious 
experience, tend to bo those cases in which one expects the tether to bo brought 
into tension, as by gravity gradient or drag farces. The reason is simple: 
SLACK3 assumes that all the segments joining it's component mosses are slack, 
except at infinitosimally short "bounces." When the physicul situation requires 
that the (continuous, physical) tether becomes taut and stretches for some 
finite time, the program can only try to approximate this by placing tho 
ball-and-spring bounces clooer and closer together. Sometimes it appears to 
make a reasonable approximation , as when the tether whips past tho Shuttle and 
rebounds at the opposite extension, but other canes simply prove too difficult, 
(One might also harbor doubts about these seemingly well handled cases, although 
the physically appealing behavior of the simulations lends some confidence.) 

Initially, we had hoped that it might be possible to handle the cases in 
which segments remain in tension for finite periods, within the same 
ball-and-spring model, It now appears that adequate numerlcal/anelytical 
handling of the tensioned case would not only be a major programming task but 
would bo computationally expensive. The number of masses required for realistic 
simulations (thirty or more) and the large number of bounces (into and out of 
tension) occupy substantial, though not prohibitive, computer time when the only 
calculations are simple root finding and free trajectory calculation. Taut 
tether calculations, which would involve either integrating sets of differential 
equations directly or elgensolutions to large (though sparse) non-symmetric 
matrices which would keep changing in form and size as segments come into and 
out of tension, is likely to be much more expensive. Some Imprecise, but 
perhaps reasonable, approaches to including taut segments at least crudely have 



Pago 30 


boon suggested and may be tried if they seem feasible after further examination. 

It is interesting to note the following! our assumption of a completely 
free tether recoil loads to a configuration on loss of tension in which the 
tether is travelling with a uniform velocity and is axtonded to just its natural 
length. Immediately, at one and it interacts with the doploymant boom, but 
throughout the tether the gravity gradient forces act to stretch the tether and 
bring It uniformly, immediately and simultaneously into tension. (This was a 
primary influence behind introducing randomization into the simulation.) Indood, 
it should have boon tensioning the toth f immediately behind the wave of 
dotonsioning. Physically, of course, we know this will not happen in such a 
clean and immediate fashion} this is an artifact of turning on the forces only 
after the idealized elastic recoil. The initial loss of tension process 
seemingly needs further study with the gravity gradient, and possibly Coriolis 
and drag, forces taken into account. 


2.3.2 Tether Properties And implications - 

Viscous damping in the tether has not been given the attention it deserves . 
In general, thia damping determines a wavelength below which tether oscillations 
cannot be sustained. More particularly, with regard to the SLACK3 model, it 
also defines a segment size below which the damping forces dominate the 
ball-and-spring bounce forces which SLACK 3 primarily deals with; Including 
damping approximately is insufficient if it is so severe that a pair of masses 
will ‘ ever exit from a bounce between them. As we shall see below, it appears 
that for soma not untypical cases (e.g. JO segments for a 500 meter tether 
remnant) the segment length is comparable to the wavelength for critical 
damping. 
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To have a handy intuitive summary figure , one related to the internal 

dynamics of the tether which are the current focus of our work, define a 
"wavelength for critical damping" of longitudinal tethor oscillation modest 

a 7 t Cv / s/ /iAE , This is found by applying the damped wave equation to a 
system with fixed boundary conditions at 0 and L, performing the standard 

separation of variables, and examining the behavior of each mode. The modes are 

sinusoids having wavelengths W n a 2L/n; those modes for which W n > Wcm are 
oscillatory, though damped! those modes for which W„ < Wcrit ere 

non-oscillatory , showing only the exponentially damped behavior. Note that this 
critical wavelength depends only on the tether properties , not on the length of 
the system. Another way of thinking of this critical length is that a freely 
hanging tethor of length less than about L/2 (wo have not actually dono thi.‘ 
calculation yet) should possess no oscillatory longitudinal modes. Attempts to 
include this internal damping (as opposed to external damping such as motion 
through a fluid) in the calculations for lateral ("plucked string") oscillations 
have resulted in a term coupling with the longitudinal oscillations. The 
effects of this term are not clear and are still under study, 

We have two sources of information as to the damping in the tethor 
material, both provided by Martin Marietta Corporation. These sources appear to 
differ substantially. First, there is a Juno 1983 report ("Analysis of Tethered 
Satellite Orbital Dynamics for Selected Mission Profiles, by Bodley and Park) 
which uses a value Cy = 9.19 x 10 4 at one point in its analysis (p, 54); the 
source of this number is not clear, and it may be a preliminary value. Second, 
an internal report on "Tether Testing' 1 dated September 1983 was kindly provided 
by A. C, Park of MMC; this report gives a "per cent damping" for an 

experimental setup consisting of a mass hanging on the end of a length of 
tether. We have calculated the damping coefficient from the experimental report 
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(bat see the caveats ■ ow) , and converted the earlier numbor to tho same (mkn) 
units, In summary; 

- Prom the experimental data, AE # !0 5 kg-m/s 2 and C v ^ 200 kg-m/s 
("averaging" the experiments described) , These imply a critically damped 
wavelength for tether longitudinal oscillations of about 25 meters (say, 
between 10 and 100) . 

- The value Cy = 9.19 ;; 10 4 kg-km/hr given in the Orbital Dynamics ..." 

report, when converted to mks system becomes 2.6 x 10 4 kg-m/s. This is 
about 100 times greater than the value calculated from the experimental 
data above, and implies a critically damped wavelength of about 3 km. (The 
values for AE and fx, converted to mks, are consistent with the experimental 
AE and assumed 3 mm Kevlar at 1.5 g/cc density.) 

Our interpretation of the experimental results, however, is somewhat 
uncertain since (1) the precise testing and computation procedures used were not 
clear, and (2) neither the direct experimental results nor a derived material 
property (Cy) were reported, hu.t a "par cent damping", presumably the damping 
ratio for the experimental setup consisting of a 4.5 kg mass on tho end of an 11 
or 21 m "spring." 

These calculations have been sent to A. C. Park at MMC for his comments and 
any light he may be able to shed on them. In any event, the implied critical 
damping scale, even of the smaller C v value, is still significant for our 
studies and deserves further experimental determination. ~ 
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2.3,3 High Resolution Loos-Of-Tension Model - 

There has been no effort on this task ir the present reporting period. 
Desirable extensions are the correct Inclusion of damping (debugging a feature 
already coded) , and the inclusion of gravity gradient force (which would retain 
the one-dimensional nature of the problem) and possibly Coriolis and drag, 
following the discussion of Section 2 .3, 1,3, 


2,3.4 Analytical Studies Of The Slack Tether Problem - 

There has been insufficient effort on this task in the present reporting 
period to report. Likely directions for continuing investigation Include: 
transverse motions; tether constitutive relations, in particular determination 
of realistic constitutive relations for single-component and multi-component 
tethers from experimental data in hand; calculation of characteristic loop size 
for a "collapsing" tether after a break. This latter will have direct 
significance to tether safety Issues. 


2.3,5 Concluding Remarks - 

The slack tether simulation program SLACK3 has been extended, made more 
easy to use, and made more resistant to failure. Several case studies have been 
simulated and analyzed. The relative effects of drag, Coriolis force, and the 
rebound of th* tether after passing the Shuttle are beginning to become 
apparent, Cases in which SLACKS has difficulty have been examined and found to 
stem primarily from physical situations where an actual tether would tend to 
come fully taut. These are Intrinsically difficult for SLACK3 to handle with 
its current completely slack segment model. Future work with SLACIC3 should 
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include minor debugging and running of cases of operational interest to NASA, 
Somo effort may be given to modifying SLACK3 to operate at least sufficiently in 
taut cases, though a complete solution to this problem seems prohibitively 
difficult, An analytic study of the situations in which the cut tether fails to 
fully recoil past the Shuttle due to gravity gradient and other forces is 
underway and should clarify, or at least allow one to predict, those cases 
causing trouble to SLACKS, 

The influence of tether properties, particularly the viscous damping Cy, on 
slack tether modeling has been examined. £xperimental data appear to be 
preliminary, but those few available span the range in which the damping may 
seriously affect both the cut tether's physical behavior and our ability to 
model it. Definitive experiment is highly recommended. 

Analytical studies are expected to progress in the next reporting period. 
Anticipated results are an estimate of the size of loops formed when the 
recoiling tether "buckles" in analogy with a column. The loop scale has obvious 
Implications for tether safety studies. The high-resolution ball-and-spring 
model for tether breaks, which complements the analytic studies, may bo enhanced 
to include damping and gravity gradient forces; out of line forces such as 
Coriolis, drag and bending stiffness might also be included. 
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2.4 Numerical Calculations Of The Electric Eield Around An Electrodynamlc 

Tether 

2.4.1 General - 

The evaluation of the potential hazards associated with the operation of an 
electrodynamic tether on hoard the Shuttle Orbiter requires a thorough modeling 
of the plasma environment in which the tether is embedded. As an initial step, 
we have embarked on an analytical effort under the simplifying assumption that 
the tether is deployed in vacuo. This initial effort will provide us with 
reference values against which we will compare the distribution of the electric 
field around the tether when plasma parameters will be introduced into the 
analysis , 

In the first quarterly report we calculated the electric field at the tip 
of a 20 kilometer long wire, moving through the earth's magnetic field, modeled 
as a 40 kilometer long prolate spheroid in vacuo. The assumption was made that 
the spacecraft to which the other end of the wire is attached could be modeled 
as an infinite ground plane. It was shown in that report that the spheroidal 
model (which has a radius of curvature of 0,5 A (5 x 10' 9 cm) at its ends, 
therefore a dimension less than the diameter of an atom) is not a realistic 
representation of the electric field near the end of a physical wire. In fact, 
that model predicted an electric field of 1.145 teravolts/meter at a distance of 
1 A from the end of the wire and 4,5 teravolts/meter at the end of the wire. It 
was estimated in the first quarterly report that the electric field predicted by 
the spheroidal model was high by a factor of 2 x 10 7 . This estimate was based 
on the assumption that the end of the wire was spherical and of the same 
diameter as the wire. To determine what the electric field would be, still in 
vacuo, at the end of a wire modeled as a uniform circular cylinder with a planar 
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end, wo have started to adapt a program which was originally written to 
calculate magnetic scalar potential with a closed boundary to calculate electric 
field with an open boundary. The mesh on which the program iteratively solves 
Laplace's equation was replaced with one in which the mesh spacing varies 
geometrically. Thus in the neighborhood of the wire the axial lines are spaced 
0.1 mm apart while at the outer edge of the mesh the spacing has increased to 
several meters. 

Since the prolate spheroid is an excellent model of the far field of a thin 
wire, we have chosen to initialize the grid with potential values calculated by 
the spheroidal model . 


2.4.2 Software Program Development - 

To minimize the numerical problems encountered in the modeling of a long 
wire with an aspect ratio of 2 x 10 7 we tested the new program with a 400 meter 
long cylinder of 0.33 meter radius in vacuo. The grid on which the potentials 
are calculated extends axially from a point 44.25 meters beyond the end of the 
cylinder and radially from the axis of the cylinder to a distance of 40.15 
meters from its axis. This distance resulted from the choice of a test grid 80 
mesh boxes long and 40 mesh boxes in radius; with a mesh width of 10 cm at the 
end of the cylinder and a geometric factor of 1.1. Thus if the mesh box at the 
end of the cylinder is A 0 cm on each side, the mesh size n grid lines from the 
end will be 1.1 n Ao, Naturally only mesh boxes which lie on a.- diagonal lino 
extending from the end of the cylinder at an angle of 45° to tho axis will be 
square. A list of the axial and radial coordinates of the grid and the mesh box 
lengths and widths is given lri Table 2.4.1. 



I 1‘ 


Pago 45 


'• A R-t a 


I Hi X 

JHA x ^oour, 

IT6AK ITHfN 

11 

»0 1 

1*00 70 

143 5 3 

U6E5F 

8913 1SPLC 

srcsr 

2' 3K 

.0013 6 

5 

3H?7MC 0C.7A2 

1 

159..VH3 

3.7404 

2 

If 3.5*87 

7.4004 

3 

166.7-31 

3.071 5 

4 

171.5774 

2.9112 

5 

1 1’, £176 

2.55*» 

6 

175. 4674 

2.3225 

7 

1 /7. 7740 

2. Ill* 

9 

179. fl -*6j 

1.9174 

7 

1-1.61*7 

1.74*9 

1 3 

1 '3. 5636 

1.5863 

11 

l J E . 1 367 

1.4421 

12 

l '6.6791 

t .31 10 

13 

18 7. 8? 31 

1.1918 

14 

1-7.0618 

1.09J5 

15 

1 '3.1553 

0.7850 

16 

1 1.1=33 

0.8859 

17 

1-2.0457 

1.8190 

!•* 

17?.«**7 

9.7450 

19 

19 J, 5-* 98 

0.6728 

’1 

1-19.2775 

1.5116 

21 

114.6 191 

1.5560 

22 

1 16.1*01 

0.5064 

25 

p-;.i4-« 

9.4595 

24 

1 "-6.47 51 

1.41/7 

25 

l c *f .8?. 18 

0.J77T 

36 

197. 2925 

3.34-7 

27 

1*7.54 77 

1.313° 

29 

1 17. 661 6 

0.2963 

29 

1*8. 1417 

3.2594 

30 

176.40*3 

0.235° 

U 

t <8.4420 

5.2149 

32 

1 '8.6664 

0.154* 

33 

171.0513 

0.1772 

34 

t’9.2294 

0.1611 


35 

r 3.36*9 

0.1464 

15 

2.79/8 

36 

175.5760 

0.13*1 

16 

1,17 72 

37 

155.6690 

0.1210 

17 

f.l’iO 

38 

159. 7060 

0,1101 

18 

4 . C j-tC 

30 

l-V.IPOO 

0.100? 

11 

4.5099 

40 

200.0000 

0.100? 

20 

5.1159 

41 

2CC. 101C 

o.nor 

21 

5.7275 

42 

2< 1.21 CO 

0.1210 

22 

6.9913 

43 

2C0.3M0 

0.1331 

23 

7.1413 

44 


0.1464 

21 

7, '843 

45 

2. C. Cl 15 

0.1611 

23 

2.0497 

46 

270,7716 

0.177? 

26 

9.8*17 

47 

2C0.54f-7 

0.1949 

27 

11.51-2 

48 

211.1436 

0.214* 

28 

19.1111 

49 

201.3690 

0.23=9 

29 

13.491} 

50 

2v 1.51,38 

0.2594 

30 

14.0411 

51 

2 C 1 , f 5 3 1 

1,2853 

31 

15.44 11 

52 

2:2.1364 

0.3138 

37 

18. 1744 

83 

S' 2. 4*23 

0.3452 

3! 

Cl. 11 3* 

64 

2.2.7-75 

C.J757 

31 

22.2952 

55 

79 3. 177? 

0.4177 

*1 

54,04 77 

56 

2,3. = ** 0 

0.4596 

36 

27. 1323 

57 

254.0545 

0,5054 

37 

27.9127 

58 

211. 5= 99 

0.5560 

3a 

33.00*0 

so 

2*5.1159 

0.6116 

39 

16. 4044 

60 

21 5.72 75 

0.6728 

40 

47.1448 

61 

2O6.4C02 

0,7400 

1 = 

-199, 517=0 

6? 

21 7. 1 4" 1 

0.8140 

?l* 

18.21237 

63 

2*.7.65*3 

0.6=54 



64 

2C8.8471 

0.9B50 



68 

215.SM7 

1,0935 

' 9 Mi r S 


66 

210.51,-2 

1.191° 


7 , 2 1 9ft 6 r *03 , 

67 

217.1100 

1.3110 

1 A C 1 12 3436799 

6S 

213.42)0 

1 .4421 

*f*t 


6? 

214. f 63 1 

1 .5PC- J 

( J pc 3 

.0) 

70 

216.44-*4 

1.7449 

*87PMT 


71 

218.1-43 

1 .9194 

< I 4 t 1 y « 

T*t OOF 5.2) / 

72 

2.-;. It 37 

2.1114 

Z^MAa = 

6. 7.0 J31-+02 

73 

1 

2 , 3 2 7 = 

r «MA' = 

...5l:3:-*02 

74 

? -4. 64 >f 

2.594* 

••■'1 Cl 123436759 

/5 

i-’.j* o* 

*.8113 

* FM r 


74 

'• !> 

i.j-r 

< IFF 1 

.21 

/ T 

: * . ' 

T , * . r t 




; ■ . v ’ 

’ 4 ?L 4 



7 

*• 

*.n«* 



^ j 


* , F * r 




- to 

1-.-48 



1 


1.1C0* 



3 

r.i::i 

'-.111- 



1 

G .21 10 

0.1210 



4 

6. 3* 1C 

0*1331 



5 

0.4*41 

0.1464 



6 

1 . f 1 C 5 

0,1611 



7 

0.7 no 

C .1772 



B 

-.9**7 

0.194* 



4 

1.14 36- 

0.2144 



10 

1.3579 

0.2358 



11 

1.5-- ’7 

C . 2 5 9 4 



12 

1.E5 31 

0.2853 



13 

2.13=4 

0.3138 



14 

2.4523 

J1..J45? 





Table 2 , 4*1 Coordinates and mesh box dimensions for the grid. 
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2,4.3 Interpretation Of Plots Generated By The Code Laplace - 

The grid whs initialized with the spheroidal function on the assumption 
that the electric field !,n which the cylinder is immersed is otherwise uniform 
and of magnitude 1 volt/meter. The region occupied by the cylinder from mesh 
box 1 through mesh box 40 in the e.xir-1 direction and from mesh box 1 through 
mesh box 3 in the radial direction is then set to zero to represent the zero 
potential cylinder. A plot of the electric equipotential surfaces around tne 
end of the wire is shown in Eigure 2.4.1. The grid which represents a region 
84,404 meters long of radius 40,145 meters is actually 22.85 cm long. The two 
equipotential lines at the top are spaced 2 cm apart on the axis. Since the 
equipotevitial lines are spaced 10.98 volts apart this indicates an average axial 
electrical field E z of 1.486 volt/meter. The calculated electric field for this 
region is shown in Table 2.4,2 from which we see an axial electric field of 
-1,43 volt/meter at position 1=0, J=74. The electric fields printed in the 
figures are calculated at the center of the edges of the mesh boxes as shown in 
Eigure 2,4.2. 

The maximum electric field calculated for position 1=3, J=39 is E r - -620,33 
V/m, E z = 819.05 V/m. This value is not very accurate, because the potentials 
on the cylinder were set to zero whiJ.e the potentials external to the cylinder 
were calculated in the absence of the cylinder. A more representative electric 
field is given at position 1=0, J=40 , when E r = 0 V/m and E z = -121.03 V/m. 

U si ng the calculated potentials on_ -the -top and bottom edges— of- the- grid,- — 
the outride radius and the surface of the cylinder as boundary conditions, the 
potential at all other points were calculated numerically using the Llebermann 
net procedure. A plot of the resulting electric field is shown in Eigure 2.4.3. 

The axial electric field at position 1=0, J=74 is shown in Table 2.4.3 as E* = 
1.55 V/m. The maximum electric field again occurs at position 1=3, J=39 , where 
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Table 2,4.2 Values of electric field computed on the basis of 
calculated potentials. 


Pago 40 


analytical ly 



Pago 49 


V(I,J+1) 

i~ 


} 




V(I+1,J+1) 

n 

i 

i 

i 

• 

i 

i 


vci,J) 




J 

Vfl+l.J) 


figure 2.4.2 This diagram shows the position at which E r and E r are calculated 
for each mesh box. The naming convention adopted is that each electric field is 
labeled by the coordinates (I,J) of the lower left hand corner of the mesh box. 



Pages 50 


E r ® **235 <69 V/m and E ( * -187.03 V/m, This result in eorract because the zero 
cylinder potential wan one of the boundary conditions for tho numerical 
calculation. Tho axial field at 1=0, J=40 is shown to be -115.76 V/m and at 
1=0, J=39 as E* = -124,35 V/m, 

The plots and tables of electric field shown here were calculated to test 
the new software program in a relatively easy case. When a wire modeled as a 
cylinder 4000 m long and 2 mm diameter was calculated, wo found that roundoff 
error caused numerical instability in tho analytic functions used to initialize 
the grid. These functions are correctly calculated in double precision on a 
Control Data computer (29 places) but the shorter word length of tho VAX does 
not permit their use. We are testing approximation techniques to calculate the 
function 



in the limit as £ approaches 1. 

This modification of the software program (which is not yet complete at the 
time of this writing) will be illustrated in the next quarterly report. 
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Table 2.4.3 Values 


of electric 


field computed on the basis of numerically 


calculated potentials 
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3.0 PROBLEMS ENCOUNTERED DURING REPORTING PERIOD 

Significant progress has been made in the development of the 
three-dimensional eloatrodynnnic computer code for the evaluation of the 
electric potential around a severed tether in vacuo. Numerical results for a 
tost case tether severance with a short tether length are presented in this 
quarterly report, Longer tethers are still difficult to deal with because of 
numerical accuracy problems. These problems have impeded us to show plots of 
the electric field around a severod tether of substantial tother length, 
However) a solution to this problem has boon singled out and computations for 
longer tethers will bo performed during the next reporting period. 

4.0 ACTIVITY PLANNED EOR THE NEXT REPORTING PERIOD 

The study of retrieval control algorithms will be continued with soma study 
of the stability of simple rate control laws. The study of thruster control 
algorithms will be continued with possible topics including control of 
out-of-plane librations, yaw control, use of the in-line thruster, and the 
effect of rigid body rotations of the subsatellite. 

Future work with SLACK3 will include minor debugging and running of cases 
of operational interest to NASA. The capability to accelerate the Shuttle in 
order to avoid the recoiling tother will be implemented in the computer code and 
cases will be run. Analytical studies on tether detensioning are expected to 
progress in the next reporting period. Anticipated results are an estimate of 
the size of the loops when the recoiling tether "buckles" in analogy with a 


column, 
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The computor code for electric field computations will bo improvod. The 
electric field around tethers in. vacuo severed far from the Shuttle will be 
computed. Modification of this computer code to include plasma characteristics 
appears to be much more complicated than expected. During the next reporting 
period simplified computations with the plasma will bo attempted to have 
indications on the potentially hazardous situations. 



